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' We discuss the resummation of the large logarithmic terms appearing in the heavy quark effects 

on parton distribution functions inside the virtual photon. We incorporate heavy quark mass effects 
by changing the initial condition of the leading-order DGLAP evolution equation. In a certain 
kinematical limit, we recover the logarithmic terms of the next-to-leading order heavy quark effects 
obtained in the previous work. This method enables us to resum the large logarithmic terms due 
to heavy quark mass effects on the parton distributions in the virtual photon. We numerically 
calculate parton distributions using the formulae derived in this work, and discuss the property of 
ClY the resummed heavy quark effects. 
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I. INTRODUCTION 



The Large Hadron Collider (LHC) Q has restarted at the CERN for the purpose of discovering the Higgs boson, 
the new physics beyond the standard model, and investigating the detailed information for the quark gluon plasma, 
the B meson decays and so on. The precise measurement will be needed to confirm the discovery of the Higgs boson 
t— I . and searching the beyond standard model at the electron-positron collider like International Linear Collider (ILC) @ 
and Super KEK-B 3]. In such a case, we have to know the behaviour of quantum chromodynamics (QCD) at high 
energies because of the largeness of QCD corrections. 

There is a well-known fact that the cross section of the two-photon processes e + e~ — > 7*7* — > e + e~ + hadrons 
dominates over that of the one-photon annihilation processes e + e~ — > 7* — } hadrons in the electron-positron collisions 
at high energies j4j]. Let us consider the two-photon processes where both of the outgoing e + and e~ are detected 
and one of the virtual photons is far off-shell with mass squared q 2 = —Q 2 , while the other photon is close to the 
on-shell with mass squared p 2 = —P 2 . In this kinematical region, the former photon is called the 'probe photon' and 
the latter one is called the 'target photon' (see Fig.[T]). 

We can regard this two-photon process as the deep-inelastic scattering in the electron-positron collision where the 
target is a photon rather than a nucleon. In this point of view, we can define the photon structure functions as 
the analogues of the nucleon structure functions, and the photon structure functions are predicted by the quantum 
electrodynamics (QED) and QCD. One of the difference between the nucleon structure functions and the photon 
structure functions is the target mass (P 2 ) dependence. The P 2 is not fixed for the virtual photon case, on the other 
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FIG. 1: Deep inelastic scattering for virtual photons in e + e collision. 



hand, P 2 is fixed for the nucleon case. There are manygood reviews for the theoretical as well as the experimental 
works for photon structure functions, for example, see [5|. 

The real (P 2 = 0) unpolarised photon structure functions F^ix, Q 2 ) and F2(x, Q 2 ) were investigated in the parton 
model (PM) |6j, in the perturbative QCD (pQCD) based on the operator product expansion (OPE) [7[ supplemented 
by the renormalisation group (RG) method @, Q , and also on the QCD improved PM powered by the parton 
evolution equation [liTfl4l |. The real polarised photon structure function pH(x, Q 2 ) was investigated with pQCD in [l5[ 
for the leading-order (LO), and for the next-to- leading order (NLO) in [la. Il7j. 

The virtual (P 2 ^ 0) unpolarised photon structure functions F^{x, Q 2 , P 2 ) and F^{x, Q 2 , P 2 ) were also investigated 
by [T^-22] in the kinematical region 

A 2 < P 2 < Q 2 , (1) 

where A is the fundamental QCD scale parameter. The advantage to study the virtual photon target for the kine- 
matical region ((T|) is that we can calculate the whole shape and magnitude of the photon structure functions entirely 
by the perturbative method. Based on the recent results for the three-loop calculation for the photon-quark and 
the photon-gluon splitting functions (anomalous dimensions) [3, US H3 , the unpolarised virtual photon structure 
function F^{x,Q 2 ,P 2 ) {Fl {x , Q 2 , P 2 )) was studied to the NNLO (to the NLO) [HHl], and the polarised virtual 
photon structure function gJ(x,Q 2 ,P 2 ) was studied to the NLO in |17l I27l - l29j . 

In the parton picture, the photon structure function is expressed as the convolution of the parton distribution 
function (PDF) in the virtual photon with the coefficient functions in the OPE formalism. We can also give the 
definite prediction for the PDFs inside the photon, which we call "photon PDFs" for short, in the present paper. 
The theoretical calculations were done for unpolarised and polarised photon PDFs in (20l I30l - l34j ]. However, these 
calculations in [3, US [IS [34| were assumed that all quarks in the virtual photon are massless. When the centre of 
mass energy is enough large to produce heavy quarks (with mass m), i.e. (p + q) 2 > 4m 2 , the heavy quark mass 
effects should be taken into account. In the case of the nucleon target, the heavy quark mass effects were studied by 
a method based on the OPE in [35| . 

Many authors have investigated the heavy quark mass effects in the photon structure functions [H, E3, HH US 13(1 - 
[42[. The heavy quark mass effects were included in the theoretical calculation for the PDFs of the real photon in 
Ref. [Ill, [36| by changing the initial condition of the DGLAP equation, for the virtual photon PDFs in Ref. (42l . \AS 
by using the OPE formalism supplemented by the mass-independent DGLAP equation. In Refs. j32j 
quark mass effects are incorporated by changing the initial condition. On the other hand, in Refs. [45 
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quark mass effects are included by evaluating the finite matrix elements of the heavy-quark operators between the 
photon states. Recently we have found that the DGLAP equation with modified initial condition for the heavy quark 
PDF [32|, |36[ leads to the similar results which we have obtained by the OPE method [42|, |43| as we have mentioned 
before. 

In the present paper, we adopt the alternative way to treat the heavy quark mass effects on the virtual photon 
PDFs by setting the initial condition for the heavy quark PDF in the DGLAP evolution equation, which amounts to 
sum up the large logarithmic terms. Through this prescription we try to improve the previous results for the heavy 
quark mass effects. 

In the next section, we discuss the basic formalism of the evolution equation for virtual photon PDFs. We derive 
the explicit expressions for the virtual photon PDFs with the heavy quark mass effects in which the large logarithmic 
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terms are resummed in section Hill In section HVl we present the numerical calculation of the results for photon PDFs. 
The final section is devoted to the conclusions. 



II. BASIC FORMALISM 



First let us discuss the basic formalism of the evolution equation for the PDFs in the virtual photon to the leading- 
order (LO) in QCD. We consider the system with n/ quarks and decompose it into two groups, rif — 1 light (i.e. 
massless) quarks and a heavy quark. Let 



qi(x,Q 2 ,P 2 ), ql(x,Q 2 ,P 2 ), G 7 (x,Q 2 ,P 2 ), r 7 (x,Q 2 ,P 2 ), 



(2) 



be light quark distributions (with i flavour and i = I, . . . ,n/ — 1), heavy quark distribution, gluon distribution and 
photon distribution inside the virtual photon, respectively. All the PDFs q] (x,Q 2 , P 2 ) evolve from the virtuality of 
the target photon P 2 to the virtuality of the probe photon Q 2 . At the LO of the QED coupling constant (O(a); 
a = e 2 /4ir), r 7 (x, Q 2 , P 2 ) does not evolve with respect to Q 2 and therefore we set r 7 (:r, Q 2 , P 2 ) = 6(1 — x). Since 
light quarks are distinguished from other quarks only through their electromagnetic charge, it is good to change the 
flavour basis to the light singlet and the light nonsinglet. We define the light singlet q 1 Ls and the light nonsinglet q~L ns 
by the equations 
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iLns 



Ql 



n f - 1 



(3) 



where is the electromagnetic charge for i-th flavour quark in the unit of proton charge. The photon PDFs are 
described by a row vector q 1 which satisfies the inhomogeneous DGLAP evolution equation 



79 >,Q 2 ,P 2 



dhxQ 2 

where the row vector q 1 (x, Q 2 , P 2 ) is defined as 



q(y,Q 2 ,P 2 )P(-,Q 



k(x,Q 2 ), 



- Inl 



ills' QhiG" 1 iQlns) ' 



(4) 



(5) 



and another row vector k^(x,Q 2 ) = (fc^ s , fc# , fc^, ki, ns ) denotes the photon-parton splitting functions. The 4x4 
matrix P(z, Q 2 ) is expressed as 



P(z,Q 2 ) 



fPE L (z,Q 2 ) Phl(z,Q 2 ) Pgl(z,Q 2 ) 

Plh{z,Q 2 ) Phh{z,Q 2 ) Pgh{z,Q 2 ) 

Plg(z,Q 2 ) Phg(z,Q 2 ) Pgg(z,Q 2 ) 

\ 



r>NS 
LL 



(6) 



(z,Q 2 )J 



where each element Pab means a splitting function of B parton to A parton. 

One can solve the DGLAP equation @ by introducing the moments of the photon PDFs Q. Here we discuss the 
procedure briefly. By taking the moment, one obtains the equation 



d\nQ 2 

where an n-th moment f(n) of a function f(x) is defined by 



q<{n,Q 2 ,P 2 ) - k(ji,Q 2 ) + q^(n,Q 2 ,P 2 )P(n,Q 2 



(7) 



fin) 



dxx n ~ x f{x). 



(8) 



Then we introduce the variable t as [44 1 



A in 

/3o a s (Q*y 



(9) 
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instead of Q 2 . We expand fc„ = k(n,Q 2 ), P n — P(n,Q 2 ) and q^ = g 7 (n, Q 2 ,P 2 ) in powers of the QED coupling 
constant a as well as the QCD coupling constant a s as follows 



< = <£ (0) +<tf (1) + -", (12) 

where the a s dependence appears in implicitly with the notation (2||, and qZ and correspond to the LO 

and NLO solutions, respectively. One finally obtains the LO solution q« for the DGLAP equation ((7|) as 

*•» - ^^^ih^}^ ^^' (13) 

y J i 1 i 

where r is the ratio of QCD couplings which is defined by 

«.(0) a s (F 2 )- 1 J 

Now the last term of fl X3[) is determined by the initial condition. Although we usually set qZ^°\o) — 0, as we took 
in our previous paper [43| . we change the initial condition for the heavy quark PDF and we will discuss the relation 
between this modification and resummation in the next section. 

00 and 0i appears in the perturbative expansion of the QCD running coupling constant of the QCD running 
coupling constant 

a s {Q ) = -0o — -. 01 ,. rn H , (15) 



dlnQ 2 ' ' " 4tt ' l (4tt) 2 

with 0o = 11 -2n//3 and 0i = 102 - 38n//3. The relation between and K,^ 0) is given by, 

W = Wn\ (16) 
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where the elements of row vector = (k^, K^f 1 , Kq", K^"^J are evaluated to be 



K°£ = 24(r V -l)( e 2 ) L fc°, (17) 

= 24e 2 H kl (18) 

K° ,n = 0, (19) 

= 24(71/ - l)((e 4 ) L - (e^Dfc", (20) 
_ n 2 +n + 2 

n(»+l)(n + 2)' 1 J 

where the charge factors for quadratic term (e 2 )^ and for quartic term (e 4 )_L are defined by 

n f — 1 n / — 1 

f — 1 * — ' rif — 1 z — ' 



Note that K L '™ + K^ n — KJ"", and K^' n is the usual flavour-singlet anomalous dimension for n/ massless quarks. 
The relation between P n "* and P" is given by 



Pi 0) - -k 0) =~ E (23) 



4 <» 4 
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where 7^ is the one-loop hadronic anomalous dimension matrix, A™ = 2/3od™ (i — tp, ±, Lns) are the eigenvalues of 
7^ which are given by 



0,n 

4> n ins "fijj-ip > 



\ n \ n 

A,;, — A T », o — 

1 

2 



0,n 



■ O.n 1 // O.Ti 0,n\9 1 ^ 0,n 0,n 

+ 7gg ± V (7^ " 7 GG ) + 4 7,/,g7g</, 



The P™ are the projection matrices in the spectral decomposition for 7^ which satisfy the following relations: 



pn pn _ r nn 



/jPj — 14x4, 



(24) 
(25) 

(26) 



where i, j = ip, ±, Lns runs over the eigenvalues, and the explicit expressions are given by 



i _i 
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•<H> A T 
1 0,n 
ri f 7i/iG 



0,ra 

7gV 

0,n 

7gV 

0,n 

7gg 



o 
















(27) 



(28) 



The most important part in the present paper is qZ^ \o), which is the initial condition for the DGLAP Eq. ([7]). We 
usually set this term qZ^°\o) to be vanishing at the LO. However, we will find that the resummed expression of the 
logarithm term ln(m 2 /P 2 ) which appears as the heavy quark mass effects in Eq. (65) in Ref. [43j is recovered by 
changing the initial condition of the DGLAP Eq. ([7]) for the heavy quark component as discussed in the next section. 



III. INITIAL CONDITION AND RESUMMATION 



In Ref. [43( it is shown that the additional terms for the moment of the PDF obtained by the OPE method are 
proportional to the LO renormalisation group parameters and the logarithmic terms. Therefore we can expect the 
possibility to derive the previous results by the LO evolution equation. We demonstrate it by changing the initial 
condition for the LO evolution equation. 



A. A change in the initial condition for the heavy quark 



The moment for the heavy quark component q^ n (t) can be projected after some calculations by using the Eq. (|13[) . 



4tt 



8^/?o a a (t) 



l + d»U/* K h ){ L r * ) + n f * + XI -A" V r ) 



n f n f — /n ± 



4 0), "(o)/ 



(29) 



where we denote /( A ±) = f(^+) +/(^— ) arL d m the last term should be determined later. As we mentioned 

previously, we change the initial condition of the DGLAP equation. Let us consider the following condition: 

ql iO) (t = 0) = (0,4 0) '",0,0), (30) 

where q^' n is the heavy-quark PDF evaluated at the scale t = (Q 2 — P 2 ). This modification for the initial condition 
of the LO DGLAP equation corresponds to a kind of the heavy quark threshold effect, because the evolution for the 
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heavy quark component is suppressed by this condition. The initial condition q^' n is determined by the equation 



4 0) '"(*m) = 
t-rn 



2 «.(!*) 



Setting t — t m in Eq. (|29p . we obtain the following result: 



4 0) '"(o)A 



i7T/3 



/3 a s (m 2 )' 



where the functions H n (t m ) and J n (t m ) are defined by 

4tt 



<^s (^m) 



1 1 <v°> n — \ n 



A" 

± ~ "T 



(31) 
(32) 

(33) 



(34) 



J n (j^m) 



1 d" 



E 



0,n 

7,/,,/, 



AIL - A? 



' m i 



(35) 



and here r m = a s (t m ) / a s (0) corresponds to the ratio between the QCD running coupling at the scale of the heavy 
quark mass and that of the renormalisation scale for the photon matrix element in OPE formalism. See the Appendix 
|X]for the explicit expressions of the PDFs in the virtual photon including the resummed mass effects. 

If we denote the variation terms as AqJ^ ,n (t) which are due to the heavy quark effects, then we recover the fixed 
order (NLO QCD + heavy quark mass effects) results as given by 



Aff£°>'"(t)/ 



(36) 



-2/3 Ai£, 



nj - 1 ~ a t 



AA 



n,± 



2/3 AA£, 



Aff2r w,n (*)/ 



AA^ 



AA^ (l - r d S) + Ai^ ± (l - r d ±) 
n/ - 1 



AC n , 



(37) 



AA 



n,± 



0,n 

1 7^ 



-2/3 AA£, 



An 



AC n 



n f XI- 



-2/3 AAjj, 



AC {0) - n (t)/ 



£AA 



G 



1 - r a ± 



(38) 



AA 



0,71 

7 G -0 



AJ - A T 



;r2A>AA£, 



Mo 



0, 



(39) 
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where g^'"(0)/ 8 ^ o = 2/?oAA^. The variation of operator matrix element AA^ is evaluated from the LO coupling 

ctsitm) together with the LO anomalous dimension K^ n in our method. Neglecting the finite term (large mass limit), 
the variation of the operator matrix element is evaluated as 

AAJf = --K„ In— T = -Uejj—. —7 ^hi-^. (40) 

H 2 H P 2 H n(n+l)(n + 2) P 2 v ' 

All the results in the large mass limit is consistent with that of the variation terms in Ref. [43[ except for the 
mass-independent finite terms. It is impossible to recover these terms only through the DGLAP equation and these 
difference could be considered as the scheme-dependence for photon PDFs. 

B. A certain limit 

Next we consider the expression by taking a limit (A 2 <C P 2 <C m 2 ) in order to compare our results with those in 
Ref. EsJ. 

In terms of the LO running coupling constant 

*«*> - tra&m- (41) 

the ratio r m can be written as 

q s (m 2 ) \n(P 2 /A 2 ) ln(m 2 /P 2 ) = 

a s {P 2 ) ln(™ 2 /A 2 ) ln(m 2 /A 2 ) ~ [ ' 

where e = ln(m 2 /P 2 )/ ln(m 2 /A 2 ). Considering the case of large-mass limit: ln(m 2 /P 2 ) <C ln(m 2 /A 2 ), we can set the 
region of e as e -C 1. Then we can expand H n (t m ), J n (t m ) in the modified initial condition term q^' n up to 0(e) as 

Air 

H n (t m ) = ^_^»[ e + 0(e 9 )], (43) 
Jn{t m ) = 1 + 0(e), (44) 
where we have used (1 + e) d w 1 + de. By using the fact - ; = /3o hi ^3-, we obtain the result as 

^ 0) ' n(0)/ 8^ " -h* * ln S [ 1 + ° (e) ] ' (45) 

where the order of the neglected term 0(e) corresponds to the terms like ln fe (m 2 /P 2 ), (k = 1, 2, 3, . . . ). We recover 
the results about the large logarithmic term due to the heavy quark mass effect which appears in the result by OPE 
formalism except for mass- independent finite terms. Therefore the equations (|A2I) . (|A3[) . (|A4[) . and (| A5|) which we 
derived by the modification of the initial condition for the heavy PDF in the LO DGLAP equation are more general 
forms and the large logarithm terms are resummed. The resummed terms form compact power terms (r m ) di in 

^ (0) '"(o)- 



IV. NUMERICAL RESULTS 



One can obtain the various photon PDFs from their moments by the inverse Mellin transformation. We show the 
results of numerical calculation for the heavy quark PDF q]j, the gluon PDF G 7 , the light singlet quark PDF q 1 Ls 
and the effective photon structure function F^ s . The last one, F^J S , is scheme-independent and is proportional to the 
total cross section of the two photon process. We consider the two cases which were measured in the experiments 
[45|,|46j]. The first case (i) resides in the PLUTO energy region where we regard the charm quark as the heavy quark, 
the second case (ii) is in the L3 energy region where we regard the bottom quark as the heavy quark, 



case (i) n f = 4, Q 2 = 5 GeV 2 , P 2 = 0.35 GeV 2 , m c = 1.3 GeV, 
case (ii) n f = 5, Q 2 = 120 GeV 2 , P 2 = 3.7 GeV 2 , m b = 4.2 GeV, 



(46) 
(47) 
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where A = 0.2 GeV is the QCD scale parameter. In both case (i) and (ii), we plot the parton distribution functions 
in the virtual photon with the DIS 7 scheme At\. In addition to the photon PDFs, we evaluate the effective photon 
structure function (x, Q 2 , P 2 ) defined by, 



f dxx n - 2 FZ(x,Q 2 ,P 2 ) = ^C7(n,g 2 )g7(n,g 2 ,P 2 ), 
Jo 



(48) 



F? s (x,Q 2 ,P 2 ) = F2(x,Q 2 ,P 2 ) + ^F2(x,Q 2 ,P 2 ), (49) 

where Cj(n, Q 2 )'s are the moments of the coefficient functions in OPE formalism and the index i runs over the 
related operators (quarks, gluon, photon), q](n, Q 2 , P 2 )'s are the moments of the photon PDFs, F7 L (x,Q 2 ,P 2 ) are 
the usual photon structure functions. The photon PDFs and the effective photon structure function contain the 
massless contribution and the massive contribution (mass effects). We have evaluated massless contribution up to the 
NLO in QCD, the massive contribution is considered up to the LO in QCD, and we add them together to get the 
total contributions by using the formalism adopted in this paper. 

We plot (a) the heavy quark PDF q] { (x,Q 2 , P 2 ), (b) the gluon PDF C 1 (x , Q 2 , P 2 ) , and (c) the light singlet quark 
PDF q[ s in Fig.[2]for the case (i), and the same three functions in Fig.[3]for the case (ii). The effective photon structure 
functions F^ s (x,Q 2 ,P 2 ) are plotted in Fig. |4] (a) and (b) for the case (i) and (ii), respectively. The abbreviation for 
the various predictions in the figures are as follows; 'Massless NLO' is the result without the heavy quark effect by 
using the massless OPE formalism, 'HQE' is the result with the heavy quark effect by using the usual OPE formalism, 
the 'Resum-HQE' is the result with the fully resummed heavy quark effects, and it is evaluated by the method in 
which the massless contribution at the NLO and the heavy quark effects at the LO are combined, and 'Resum-HQE- 
Exp' is the result with Eq. (|45p . respectively. The difference between 'Resum-HQE- Exp' and 'HQE' arises from the 
mass-independent finite term which one cannot reproduce only through the renormalisation group equation as we 
mentioned before. 

In general, we can see the large suppression due to the heavy quark mass effects and its resummation effect on 
the charm quark distribution (a) , the gluon distribution (b) in Fig. [5] and the effective photon structure function 
(a) in Fig. [4] On the other hand, we find a small suppression effect on the bottom quark distribution (a), the gluon 
distribution (b) in Fig.[3]and the effective photon structure function (b) in Fig.UJ For the light singlet quark PDF, the 
heavy quark effect given in Eq. Q36p is almost negligible, namely the three curves with heavy quark effects in Figs.[2Jc) 
and[3Jc) overlap with each other and they also coincide with the plot 'Massless NLO'except for the small x region in 
the case of Fig.[2jc). This is because the right-hand side of Eq. (|36|) can be written as (1 — l/n/)2 J 9oAA^-/ n (r), where 

0,n \ n O.n \ n 

/n(r) = ~(rf* + + ^3Af ( 5 °) 

which is extremely small as discussed in [43]. While the right-hand side of Eq. (|3"T)) is written as (1 — 
l/n f )2p AA™g n (r), where 

g n (r) ee (rfi + — - | l±t__ {r fl + ^t_± {r f- L = (r) «$ /(1 _ 1/n/) + /„( r )/( n/ - 1), (51) 

which is expressed approximately as (r) d ^/(l — 1/nr), and hence the heavy quark mass effects become sizable for 
heavy quark PDFs. Note that there exist no heavy quark effects on light nonsinglct quark PDF q~[ ns as seen from 
Eq. (f3"9"|) . Since each light quark PDF, q l L , is a linear combination of the singlet and nonsinglet quark PDFs, the 
heavy quark effects on the light quark PDFs are in fact negligibly small. Phenomenologically interesting features 
are the differences of the size for the heavy quark effect on the heavy quark PDF and the gluon PDF. The origin of 
this difference for the heavy quark effects comes from their electromagnetic charges. Since the absolute value of the 
charge of the charm quark is larger than that of the bottom quark, we obtain the larger reduction of the gluon PDF 
in the case (i) than that of in the case (ii). We can also see these differences in the physical observable; the effective 
structure function in Figs. |4] (a)-(b) in both case (i) and (ii). 

The resummation effect of the large logarithmic terms due to the heavy quark mass is also larger in the case (i) for 
the three functions (charm, gluon, effective structure function) than those in the case (ii). In addition to this feature, 
we observe a little bit interesting property of the (a) in Fig. @] It seems to be that the curve with fully resummation 
effect (Resum-HQE) slightly close to the experimental data than the curve without the resummation effect (HQE, 
Resum-HQE- Exp) . It might be due to the validity of our present method for the resummation of the heavy quark 
mass effects. However we cannot say anything about the validity for L3 case (the figure (b) in Fig. @| due to the small 
bottom's mass effects on the effective structure function. 
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n,=4, 2 =5 GeV 2 , P 2 =0.35 GeV 2 , A=0.2 GeV, m c =1 .3 GeV 



n f =4, Q =5 GeV , P 2 =0.35 GeV , A=0.2 GeV, m c =1 .3 GeV 



n,=4, 0=5 GeV , P 2 =0.35 GeV , A=0.2 GeV, m c =1.3 GeV 



Masslcss NLO - 
HOE - 
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(c) 



FIG. 2: Parton distributions in the photon in DIS 7 scheme for Uf = 4, Q 2 = 5 GeV 2 , P 2 = 0.35 GeV 2 with m c = 1.3 GeV and 
A = 0.2 GeV: (a) xc 1 (x, Q 2 , P 2 )|dis 7 ; (b) xG 1 {x, Q 2 , P 2 )dis 7 , and the light singlet quark distribution (c) xql s (x, Q 2 , P 2 )dis 7 • 



n.=5, =120 GeV^, P^=3.7 GeV 2 , A=0.2 GeV, m b =4.2 GeV 



n<=5, 0^=120 GeV 2 , P 2 =3.7 GeV 2 , A=0.2 GeV, m h =4.2 GeV 
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FIG. 3: Parton distributions in the photon in DIS 7 scheme for n f = 5, Q 2 = 120 GeV 2 , P 2 = 3.7 GeV 2 with m b = 4.2 GeV and 
A = 0.2 GeV: (a) xb J (x, Q 2 , P 2 ) | D is T ; (b) xG J (x, Q 2 , P 2 )dis t , and the light singlet quark distribution (c) xq[ a {x, Q 2 , P 2 )dis 7 • 



Thus, we conclude that the large suppression effect exists in the theoretical prediction for the charm PDF, gluon 
PDF in the (b) with DIS 7 scheme, the effective photon structure function (a) of Fig. |4] at the PLUTO 's kinematical 
point not because of the NLO QCD corrections, but because of the heavy quark (charm quark) effects. On the other 
hand, we can see stable results of the theoretical prediction for the bottom PDF, the gluon PDF in the (b) with DIS 7 
scheme, the effective photon structure function (b) of Fig. |4] at the L3's kinematical point despite the NLO QCD 
corrections and heavy quark (bottom quark) effects within this formalism. 
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FIG. 4: Effective structure functions for PLUTO (a); n f = 4, Q 2 = 5 GeV 2 , P 2 = 0.35 GeV 2 with m c = 1.3 GeV and 
A = 0.2 GeV: L3 (b) n f = 5, Q 2 = 120 GeV 2 , P 2 = 3.7 GeV 2 with m b = 4.2 GeV and A = 0.2 GeV. 
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V. CONCLUSION 



We have discussed the resummation of the heavy quark mass effects on the photon PDFs by changing the initial 
condition for the LO DGLAP equation. Our method is based on the mass-independent DGLAP evolution equation 
and a change of the initial condition for the heavy quark component in the LO solution of the DGLAP evolution 
equation. 

We recovered the previous results based on the NLO OPE formalism [43| except for finite terms by taking a certain 
limit: ln(m 2 /P 2 ) < ln(m 2 /A 2 ). By using the method, we can resum the large logarithmic terms due to heavy quark 
mass effects on photon PDFs in the virtual photon. We evaluate the size of the resummation effects numerically. 
The resummation effects for the charm quark PDF, gluon PDF are larger than that of the bottom quark PDF due 
to the heavy quark electromagnetic charge. We also see that the resummation of the heavy quark mass effects on 
the effective structure function tends to close to the experimental data in PLUTO case. Although we only presented 
the PDFs in the DIS 7 scheme, the PDFs in the MS scheme with our resummed heavy quark effects could be inferred 
from the difference between the PDFs in MS and those in DIS 7 in our previous paper [43[ where both schemes were 
explicitly presented. 

Now some comments on the future extensions of the present work are in order. One of the possible extensions 
is the NLO analysis. We can solve the NLO DGLAP equation including the change of the initial condition for the 
LO DGLAP equation. Other possibility is the extension of this idea to the rif quarks system with two or more 
heavy quarks. In the case of two heavy flavours, this can be achieved by decomposing the nj quarks into the nf ~ 2 
light quarks and two heavy quarks. Such an extension will be useful to analyse the virtual photon structure functions 
under the situation which contains both the massive bottom quark and the charm quark at Super KEK-B Q . Another 
example for the phenomenological application of this idea is to analyse the system with the heavy supcrpartners in 
the supersymmetric QCD. 
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Appendix A: Summary of explicit expressions for the LO solution 



We can obtain the expression with the heavy quark effects for Ls, H, G, Lns components. Each elements of the LO 
solution is defined by 



l (0) = (<& ),B (t), 4 0) '"(')> G^(t), q ^;(tj) 



(Al) 



Then the LO solutions are summarised as 
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G (0) - ,l (t)/ 



4tt 



8^/?o a a (t) 



0,ri 

7gV 



0,ra 

7 G ^ 



a 



!tt/V 



(A4) 



47T 



8tt/3o a a (t)l + d| 



(A5) 



The first term of the above results corresponds to the massless result and the second term corresponds to the variation 
term due the resummed heavy quark mass effects. Note that there is no extra contribution due to the modification 
of the initial condition for the Lns component. This is consistent with the result in our previous work [431 ] . These 
results are the new and main results in this paper. 
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